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^ ' Abstract 

' Links of factorization theory, supersymmetry and Darboux transforma- 

I tions as isospectral deformations of difference operators are considered in the 

^>«^ . context of soliton theory. The dressing chain equations for factorizing oper- 

' ators of a spectral problem are derived. The chain equations itselves yield 

nonlinear systems which closure generates solutions of the equations as well 
■ as of the nonlinear system if both operators of the correspondent Hirota bilin- 

' earization are covariant with respect to Darboux transformation which hence 

defines a symmetry of the nonlinear system as well as of these closed chains. 
I ' Examples of Ifirota and Nahm equations are specified. 

> 



1 Introduction. 

The growing interest to discrete models appeal to necessity to widen classes of 
symmetry structures of the correspondent nonlinear problems ffl. Such theories 
^ as conformal field result in set of equations to be investigated Next kind of 

' problems contain solvable lattice models @, 

Very recently a good basis for new searches in this field of differential-difference 
and difference-difference equations was discovered § in the context of classical Dar- 
boux Transformations (DT) theory development. Likely the differential operator 
case it would have links to Hirota bilinearization method and to factorization 
theory [Q with similar applications possibilities. 

One of very promising approach to the solutions construction is based on the 
dressing chain notion (lO), [0, The approach cover soliton, rational, 

finite-gap and other new solutions inside universal scheme, reducing the problem 
to a solution of closed sets of nonlinear ordinary equations with bi-Hamiltonian 
structure jist . 

In this paper we reformulate Darboux covariance theorem from introducing 
a kind of the difference Bell polynomials. Such polynomials have natural correspon- 
dence to the differential (generalized) Bell polynomials in its nonabelian version Q , 
B. Its usage shortens the transformation formulas and helps to apply the theory 



in comphcated cases of joint covariance of U-V pairs g], |16| on a way of the 
dressing chain equations derivation. 

The example of Hirota equation p^ , connected with known applications, is 
studied in Here we derive the dressing chain equations and study simplest 
solutions of it. In the last section we propose the lattice Lax pair for Nahm equations 
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p^ , see also |T^, on physical aspects. The Lax pair is covariant with respect to 
a combined Darboux-gauge transformations that generates the dressing structure 
for the equations. 

2 Darboux transformations in associative ring 
with automorphism. 

In this section we reformulate, and analyze the results from |^] in a context of its 
further use in the chain equations derivation and joint covariance of operators pairs 
investigation. We would begin with general notations. Let R be an associative 
ring with an automorphism, implying that there exist a linear infertile map T, R— >R 
such that for any ifj (x, t) , (f (x, t) gR, x e R''\t e R, 

T{^^)=T{^)T{^), T{1) = 1. (1) 

The automorphism with the defining property (Q) allows to write down a wide 
class of functional-differential-difference and difference-difference equations [|| start- 
ing from 

N 

i^t{x,t)= U.nT"'^. (2) 

For example, the operators T could be chosen as, 

Ti;{x,t) ^ij{qx + d,t), (3) 

where q G GL(n,C),(5 G i?". 

(x) = Wip{x)W-^ , W e GL{n, C) 

We would save the notations and conditions of the mentioned paper ^ discussing 
other potentials definitions in the end of the section 4. 

Let us consider the definitions of two Darboux transformations for solutions of 

D^f = f- a±r±V, 'y^^v (T±V)"' • (4) 

where if is a particular solution of the same equation (|^) For the case of a differential 
ring and for Tf (cc, t) = f {x + 6,t) ;x,S G R the limit df = lim^^o - 1)/ (x, t) 
give the link to the classical DT. 

To derive the DT of potentials Um it is necessary to evaluate the derivative of 
the elements with respect to the variable t (say - time) . We shall do it by means 
of the special functions (differential Bell polynomials analogue) introduction similar 
to [Q . Let us start from the first version of the DT definition expressing the 
Tip from (^ , hence T(p = (cr^) ^ (p. Acting to the result by T, taking into account 
(|) yields TV = T(^{a+y^^TLp = T (^{a+)~^^ {a+y^ p. Repeating the action, 
one arrives at 

m— 1 

T"<^ = n {T" (a+))-V - B+ (a+) ^. (5) 

k=0 
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Here and below the product is ordered by the index (fc) from right to left. 
Definition The equation (||) defines the function 

m — 1 
fc=0 

It is convenient to write down the t-derivative of the a by means of the functions 
5+ (ct+ ) that are connected generalized with Bell polynomials . 

N 

^t^ E [U^Bt,{a+)a+-a+T{Ura)B+^^,{a+)a+] (6) 

m=-M 

The resulting equation is some nonlinear equation associated with the 
that is reduced to a generalized Miura transformation in the stationary case (see 
sec. 3). 

The Matveev theorem gives powerful generalizations of the conventional Dar- 
boux theorem, proved originally for the second order differential equation (for gen- 
eralizations look jl^ as well) and can be formulated by means of the introduced 
entries in the following way. 

Theorem 1 Let the functions f and ip (zR satisfy the equation (|^), then the 
function — D'^ip satisfies the equation ipf {x,t) = X]m=-j\/ U^T"^ip'^, where 
the coefficients are evaluated from the following recurrent relations 

U+,, = U^M (7) 
C/+ - U+(T+ = Ui- a+TUo - at (8) 
U+ - f/^-i^™- V+ = - a+T[/„_i (9) 

U+ =a+iTUN){T^a+)-' (10) 

The expressions (|^) — (|o|) define the Darboux transformations of the coefficients of 
the differential equationl^ (potentials ) as the recurrence. The straight corollary of 
this statement solving the recurrence (^ by means of (|^) yields 

m+M 

K - E '^-^^+' - ''^ (TU^M+i-i) Btu+i (B+ (^+))"' (11) 

1=0 



U+ ^a+iTUr,){T^<J+) (12) 

The theorem establishes the covariance (form-invariance) of the equation (^ with 
respect to the Darboux transformation. 

Proof. For the proof it is necessary to check the additional equality that appears 
by the term T"^^ with essential use of the expression for from (^). ■ 

The formalism for the second DT from (^ may be constructed similarly on the 
base of the identity 

T'"(^ = n i'^') V = B-, (a-) T-\. (13) 

fc=0 
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The definition of the second type lattice Bell polynomials B^^ {a ) may be extracted 
from (p^). The evolution equation for a~ is similar to 



N 



= J2 [U,nB-^{a-)-a-T-'{U^)B-^_,{<j-)] (14) 

m=-M 

and gives the second generalized Miura map for a stationary solutions of the (0) . 



Exphcit formulas for U' are similar to (|), (|), Q, Q, (|T2|). 

3 Stationary equations as eigenvalue problems and 
chains 



The stationary equation that corresponds to ^ appears when the solutions of the 
following constraint equations ipt {x,t) ~ Xip, ipt {x,t) — ipX, are considered, or, in 
the first case, 

N 

^ C/,„T"V = AV-. (15) 

m=-M 

The derivative 

(^t = 'Pt {T(py^ - (p{Tipy^{Tip)t {Tipy^ = ^cr - a^i. (16) 

It is zero if the a and fi commute. Recall that ip G {ipx^fj,}. 

One arrives at the algebraic connection between potentials and a that is the 
corollary of either (^ or ( p^ ) in this case of (|l^). When the reduction is such as 
all the coefhcients in (|^) are functions of some unique potential u, the connection, 
allows in principle to express the potential u as the function of a. This connection 
has the same form as for the potential u as for a "dressed" one. If one takes 
connections for both ones with the correspondent sigmas and plug the result into 
the Darboux transformation formulas, the chain equations yield. 

Let us illustrate the algorithm by an 

Example of generalized ZS problem. 

Let us take the equation ^ in the minimal variant 

^t{x,t)^ {Uo + UiT)^p. (17) 

The covariance of this equation with respect to the DT of first kind means the 
invariance of Uq- We rewrite the transform of C/i by means of ( p^ 

U+ = a+ {TU,) {Ta+y' (18) 

For the spectral problem 

(Uo + UiT) if ^ flip (19) 

the connection of the unique potential Ui and is obtained form (^ and ( p^ or, 
directly from the equation (p^) , 



Ui = (/i - Uo)<J+ (20) 
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Introducing the number of iteration n for Ui, and n+1 for yields 



f/i [n+l] = {fi- Uo)<J+ [n + l] = a+ [n] T{ip- Uo)<J+ [n]) {Ta+ [n]) ' . (21) 

It could be the chain equation for the generalized ZS problem. We rewrite 
the chain equation (^) in more compact form changing the notations as follows 
Uo J, t/i — > U, cr"*" [n] an and supposing that T{J) = J. We arrive at 



(7„+i ^ J)- V„ {^l ~ J) (22) 

This dressing however is obviously almost trivial. Such phenomenon is well-known 
in the differential operator case [0. The alternative and effective transformations 
appear if in the stationary equation ( p^ one introduces the element ^ that does 
not commute with a and changes the order in the r.h.s. between the elements fj, 
and (p. (It is the second case mentioned at the very beginning of this section.) 

(J + UT)ip^ ipfi. (23) 

The formula for the potential is changed 

U = ipfiiTip)-'^ - Ja (24) 

and the role of can play the function s = ipfi{Tip)^^ . The equation (^ connects 
U and a 

Jcr + U - aJa + crT{U){Ta)-'^ = [fi,a]. (25) 

The algorithm of the explicit derivation of the chain equations begins from the 
equation ((p5|)) solving with respect to U in appropriate way. For matrix rings, it 
may be a system of equations for matrix elements, that could be effective in low 
matrix dimensions of ((p5|)), as in Otherwise it opens the special problem. Let 
us rewrite (^5|) and the DT ^ in terms of s, excluding U from denoting the 
number of iterations by index 

U[n\ s„ - J(T„. 

The equation ( p5| ) transforms as 

s - crJa + aT{s){Ta)-^ + <tJ = [/i, a]. (26) 
The use of this result give for the DT 

Sn+l - S„ = Jo-„+l + CT„ J(T„ - (T„]. (27) 

Then, taking the result (p7|)for two indices (e.g. it could be n,n-|-l) one have 
the chain system. In the section 5 we will give the explicit example for the bilinear 
Hirota equation. 

Let us mention that the chain equations for the classical ZS problem and two 
types of DT transformation were introduced in jl2|. The closure of the chain 
equations specify classes of solutions. The chain equations for a standard Sturm- 
Liouville problem were derived in connection with quantum mechanical problems 
in early Schrodinger paper (look the review [^|). In connection with the celebrated 
scalar KdV important possibility is studied in . A periodic closure of the chains 
produces integrable bihamiltonian finite-dimensional systems and, in some special 
cases, the finite-gap potentials H. 
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4 Joint covariance of equations and nonlinear prob- 
lems 

Let us consider a pair equations of the same type (j^) for a function i/j 

N 

i>t{x,t)^ J2 U„^T"'i>, (28) 



m=-M 



N' 

^y{x,t)= K„T"V. (29) 

m=-M' 

The compatibihty condition for them is the nonhnear equation 

Usy -Vst^Y. ^fe^'' - Us-kT'-'' (Vk) (30) 

for s = ~M~M',...,N + N', k € {k' ^ -M',..., iV'}n{s-fc = -Af ,...,iV }. 
In the simplest case of ZS operators in both (|2|) , (|29| ) with the subclass of 
stationary in y solutions one have three conditions 

Uot = VoUo-UoVo, (31) 
Uu = VoUi-UoVi + ViTiUo)-UiTiVo) (32) 

V^T (U,) ^ UiT (Ui) (33) 

The connection with polynomials of a differential operator and hence with classic 
Bell polynomials theory may be found if one change the definition of potentials. It 
is clear that if the automorphism T is the shift operator Tf{x) = f{x + i5), the 
coefficients of the polynomials in T should be arranged as follows. 

N in , s 

E EL™, (-ir"'^^>- (34) 

in=~M r=Q ^ ^ 

The Newton binomial formula have been used. The recursion equation that define 
classic differential Bell polynomials jsj 



r=0 



— E ( )^m-r2/r+l, (35) 



together with the definition (^ of 5+ connects the expression for these special 
functions. Let us remark that the transformations for Um found in the sec. 2 give 
the transforms for Um defined by (^). The possibility of inverse transition depends 
on the independence of functions (T — 1)"/ for a given T and the set of functions 
ip under consideration. The joint covariance of the system (|2^), ( p9| ) hence may be 
investigated along the guidelines of ||] , ||l^ , where so-called Binary Bell polynomials 
are used to form a convenient basis. 



6 



5 Nonabelian Hirota system 

5.1 Derivation of the chain equation 

Let us consider a pair of Zakharov-Shabat type equations, 

iPtix,y,t)= (Fo + l/iT)V, (36) 



i;y{x,y,t)^ (Uo + U-iT-^)^. (37) 

It differs from one used in previous section by the change T in the r.h.s. 

In a t-lattice version of the equation (0) with j S Z, one goes to 



N 

f{x,j + l)= Yl UmT^f. (38) 

The lattice in all variables case is generated by the transition to the discrete 
variables x,y,t ^ n, j, r € Z; f {x, y, t) — > /„ {j, r) , defined as in The operator 
T acts as the shift of n: T/„ {j, r) = fn+i (j, '') ■ The corresponding equations of the 
system (|6|) , (|3^) looks as 

fn (j - 1, r) = f„+i {j, r ) + v{n, j, r)fn {j, r ) (39) 



fn (j, r - 1) = /„ (j, r ) + u{n,j, r)/„_i (j, r ) (40) 

with a specified potentials. The compatibility condition of the linear equations 
^ , (|4^) has the form 

u{n,j -l,r) - u{n + 1, j,r) = v{n,j,r- 1) - v{n,j,r), 



v{nj,r - l)u{nj,r) = u{nj - l,r)v{n - l,j,r). (41) 
The second equation of (^) is automatically valid if 

u(n, j, r) = r„+i {j, r - 1 ) r^^ (j, r - 1) t„_i (j, r ) r'^ (j, r ) 

i;(n, j, r) = r„+i (j - 1, r ) t,7^ (j - 1, r) t„ (j, r ) r^^^^ (j, r ) . (42) 

We stress that the form of these expressions is valid for the order of entries that had 
been shown in (|4l|), (^2|): the form differs (generalize) from The substitution 
of the (^) to the first equation from ( |4l| ) leads to the generalized Hirota bilinear 
equation 1 1^ , compare also with generalizations in |^] 



Tn+i (j - 1, r - 1 ) r„ 1 (j - 1, r - 1) r„_i (j - 1, r ) t„ ^ {j - 1, r ) - 

r„+i (j - 1, r - 1 ))r^^ {j - 1, r - 1) t„_i (j, r - 1 ) t^^ (j, r - 1 ) - 
'rn+2 (i, r - 1 ) T^l^ {j, r-l) t„ (j, r ) t^_/:i (j, r ) + 

T„+i (j - 1, r ))T,7i (j - 1, r) T„ (j, r ) t-^, {j, r) = (43) 
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In the scalar case the system reduces to the Hirota bihnear equation ||^ 

T„ (j + 1, r ) T„ (j, r + 1 ) - r„ (j, r ) t„ {j + 1, r +1) + r„+i (j + 1, r ) r„_i (j, r + 1 ) = 0. (44) 

Using substitutions ( p^ ) and all DT formalism we might give in a nonabelian 
form that is important in connection with the applications in the theory of a quan- 
tum transfer matrices for fusion rules and the theory of quantum correlation 
functions jigt] . 

Let us return back to the DT theory. The equations (|3^) , js^) are jointly co- 
variant. Hence the theory of solution of the systems (^ ) or ([f3|) is based on the 
symmetry of the systems that is generated by the joint covariance of the equations 
(P9) , (Efl) with respect to the one of transformations (0) , namely to the 



^" (j, r) = V - a-T- V, = ^ (T- V) ' ■ (45) 



The form of the both linear equations ( |39| ) , (^0|) is obvious reductions of the 
equations ( |3^ , ( |37| ) with Vi = 1, Vq = v; Uq = 1, U-i = u. We show further some 
details in the proof of the covariance theorem because it shows the important fea- 
tures in the procedure of the chain equation derivation. Let us start, say, from (BG) . 
The conditions of the covariance are obtained as the coefficients by "0, T~^il), T ijj. 
The first one is valid automatically. 



u = u — a (r — 1) + (7 (r) (46) 



u-T-^ff- (r) ^ a- {r -l)u (47) 

Excluding the transformed potential u^one arrives to the equation that links the 
potential and the function (r) 



u T-^a- (r) - a- {r-l)u= (ct" (r - 1) - (J- (r)) T-^a^ (r) . (48) 

Note that the expressions in the equalities (|5^), (^), ( ^ are still general (non- 
abelian) and may be used in simplest (but definitely rich) closures (e.g. Cj^+i = '^n )- 
The generic is outlined in the end of the sec. 3. For some Lie algebraic approach 
to the nonabelian chains look also in iflq. 



In the scalar (abelian) case one can easily solve the equation (48) with respect 
to the potential 

_ (g- (r-l)-a- (r))r-V- (r) 

T-ia-(r)-a-(r-l) ' ^ ' 

Supplying the entries of the equation (}49|) by the index N (iteration number) 
and substituting into the equations ( ^5|) , ([f7|) we obtain two equivalent forms of the 
chain equations. For example, 

UN+i = ujv - o-^ (f - 1) + crj, (r) 

yields 
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{<^N+1 - 1) - ('')) T Vjv+i (r) _ (cr^ - 1) - CTn (^)) f^TV - 1) 



T ^<^N+i (i') - <^N+i (r-l) T ifT^Y (r) -(T^{r- 1) 

(50) 

The chain equation ( |50| ) generates the chain equation for the specific case of 
the system ( |4l| ) by the choice x n,T/„ (j,r ) = (j,r ) . A solution of the 

resuhing chain equation generates the solution of the system by use of the 
connection formula ( p9| ) and the corresponding formula for v. The transition to the 
Tn{3,r ) - functions is made by (E2h . 



5.2 On solution of the chain equation 

Let us denote 



(ct^ (r-l)-cr^ (r)) 



then the dressing chain equation (p3) reads 



Iterating this recurrence yields 



^^+1 = ^AT ^''^-" (51) 



SN+q = STV-pfg „„i _ T^- (52) 



In analogy with continuous case let us consider the periodic closure of the chain 
(|5C|), starting from the simplest case q=0, that means ctat+i = = a. As sat+i = 
Sat, 

T(T~ (r — 1) = cr~ (r) . 

It means cr~ (r + p) = T^a^ (r) . If a boundary condition in the point r=0 is 
given, then 

a- (p) = T^a- (0) . (53) 
The equation for (p{p,x) then looks 

Tip{p, x)^T {a{p, x)) <fi{p, x). 
The solution depends on the choice of T; e.g. if Tip{p, x) — ip{p, x + S), then 

ip{p, x) = exp (Ax) 

with A from 

AS ^ In [TP+^a- (0)] 
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6 Nahm equation 



Considering the following example we change a bit the DT formulas, showing the 
alternative version, similar to We would stress, however, that the formulas from 
the first section give the equivalent result. Some generalization we need within the 
reduction constraints - we need' additional (gauge) transformation denoted by g. 
This is expressed by the following 

Theorem 2 The equation 

tpy — uTtp + vip + wT~^ip (54) 
is covariant with respect to DT 

m ^ 9iT - <j)^^-\ (55) 



Where a — {T(j))l(j), (p is a solution of the same equation ( |5^ J and g is an 
invertihle element of the ring to be defined later. The transforms of the equation 
coefficients are 

u[l]^gT{u)[T{g)]-\ (56) 
v[l] = gT{v))g-' - gang-' + gT{u)T{a)g-^ + gyg'\ (57) 

w[l]^gaw[T-\ga)]-\ (58) 



Proof. The substitution of (q5|) into the transformed equation (|5j) give four 
equations assuming T"^ are independent. Three of them yield transformed poten- 
tials |5^. The fourth equation after use of the transforms takes the form 

(Ty^aF - {TF)a (59) 

where 

F ^ ua + V + w[T-\a)]-\ 

one can check the condition by direct substitution of the a definition and use of the 
equation for (p. ■ 

Remark The theorem 2 is naturally valid for the spectral problem 

AtA = uTi; + v'ip + wT~^il> (60) 

with the only correction: the last term for the transform v[l] is absent. The equation 
goes to the "Riccati equation" analog for the function a 

^^ua + v + w[T-'^{(j)Y^, (61) 

Note that the plugging of the element a = (T(j))/(j) into the equation (^ij) transforms 
it to the spectral problem for (60) (with the spectral parameter /i). 

The Nahm equations could be written in the Lax representation using the spec- 
tral equation ( |60| ) and the evolution equation 

^Py = iq+pT)^P, (62) 

The covariance of this equation with respect to the DT ((p5|)) may be established 
similar to the theorem 1 with account of the evolution for the function (j{y) 

(Ty = T{q)a — apa + T{p)T{a)a — aq = 0, (63) 

that proves the following transformation formulas for the coefficients in ((|62|)). 

p[l]^gT{p)[Tigr\ (64) 
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q[l] = g[T{q) -ap + T{p)T{<j)]g-' + gyg-\ (65) 
The joint covariance principle 1^ define the connection of potentials p, q with u, v 

p = u + pi, q^v/2. (66) 
Hence the DT-covariance means integrability of the compatibility condition of 



Uy = ^{uT{v) - vu) + (3{T{v) ~ v), 



Wy — —vw — wT — 1(f). 



Vy = uT{w) — wT u + j3{T{w) — uj), 
1 

One more possible specification is the use of periodic potentials in the problem 
^2|) with the evolution (jsj with account of the connections (|6|) that results in 
commutators in r,h,s of the equations. Some linear transformations and rescaling 

■u = a(-i(pi/2 - (ps), (67) 
V = (f3, (68) 

w^a-\-iipi/2 + ip3), (69) 

produce Nahm equations (for the periodic functions Tipi — ipi, that by the way 
does not mean a periodicity of solutions of the Lax pair ip, cj) and the corresponding 
a = {T<t>)/<P) 

ifiiy = leikiiipk,^!]- (70) 

The numbers a, (3 are free parameters. This system is covariant with respect to 
the combined DT-gauge transformations if the gauge transformation g = exp G is 
chosen as follows. 

Gy = a[i^3 + ifi/2)T{a) - a((^3 + ¥'i/2)]/2, (71) 

Finally, the following theorem may be formulated. 

Theorem The system (^) is invariant with respect to the transformations 

^i[l] = g[(^i/2 - i^3)T{g)-^ + a(^i/2 + i^3)[T-\9a)]-\ (72) 
ip2[^] = 9^2 + a[i<7ipi/2 - npiT{(j)/2 + aLp^, ~ T((p3cr))]g"\ (73) 
V3[l] = g[{^^Vl|2 - V3)T{g)-^ + a{-iv^/2 + ^3)[T-\g<j)]-^] (74) 



with the function g — expG, where G is obtained integrating (|7l|). 

Remark Similar statement may be formulated for the discrete version (see pT[) 
of the Nahm system (^0|) as it may be seen from the previous section. 

Conclusion The productivity of the method is obviously dependent on the 
possibilities to solve chain equations. The direct applications of DT formulas could 
be effective as it is well-known p7|. The approaches to the direct solution of the 
chain equations is formulated in ]13| | , and refs therein. 
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